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ABSTRACT 

We reveal the role of clustering of sub-clumps, which is expected in the cold dark matter (CDM) 
universe, in forming a bright elliptical galaxy (BEG) from a low-spin seed galaxy. This can be done by 
comparing the evolution of a low-spin seed galaxy including small-scale density fluctuations expected 
in the CDM universe (Model 1) with that of a completely uniform one (Model 2), using numerical 
experiments. We show that Model 2 cannot reproduce the properties of BEGs and forms a disk which 
is too compact and too bright due to the conservation of the initial-small angular momentum. In 
Model 1 clustering of the sub-clumps caused by initial small-scale density fluctuations leads to angular 
momentum transfer from the baryon component to the dark matter and consequently a nearly spherical 
system supported by random motions is formed. Moreover the collisionlcss property of the stars formed 
in the sub-clumps prevents the dissipative contraction of the system, leading to a large measured half- 
light radius. As a result, the end-product is quite well reproduces the observed properties of BEGs, such 
as the de Vaucouleurs light-profile, typical color and metallicity gradients, the large half-light radius, the 
small ratio of the rotational velocity to the velocity dispersion (V/a). We conclude that the clustering of 
sub-clumps, i.e., the hierarchical clustering, plays a crucial role in the formation of BEGs from a low-spin 
seed galaxy. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: stellar 
content — galaxies: structure 



1. INTRODUCTION 

Observed bright elliptical galaxies (BEGs) have the 
common interesting properties, which are expected to pro- 
vide strong constraints on their formation scenario. For 
example, Davies et al. (1983) found that the ratios of the 
rotation velocity to the velocity dispersion, V/a, of BEGs 
are much smaller than those of faint ellipticals. In pho- 
tometric properties, the light profile of BEGs is well de- 
scribed by the de Vaucouleurs law, whereas faint ellipticals 
usually have an exponential light profile (e.g., Caon, Ca- 
paccioli, & D'Onofrio 1993). 

On the other hand, Kawata (1999) showed that in the 
CDM universe, i.e., the hierarchical clustering scenario, el- 
liptical galaxies can be formed from a low-spin overdense 
region, using numerical simulation which includes almost 
all the important physical processes in galaxy formation 
self-consistently. He adopted a semi-cosmological model 
proposed by Katz & Gunn (1991). In this model, the pro- 
cess of galaxy formation is mimicked by a collapse of a 
top-hat over-dense sphere (a seed galaxy) which initially 
follows a Hubble flow expansion and has a solid-body ro- 
tation as an effect of the external tidal field. The amount 
of the initial solid body rotation is specified by a dimen- 
sionless spin parameter, A (Peebles 1971). In addition, 
small-scale density fluctuations expected in the CDM uni- 
verse are superposed on this sphere which consists of dark 
and baryonic matter. Using this model, Katz & Gunn 
(1991), Katz (1992), and Steinmetz & Miiller (1994, 1995) 
studied only the case of A = 0.08 and seemed to suc- 
ceed in forming a system which has similar properties 
to an observed disk galaxy in several respects (see also 
Koda, Sofue, & Wada 2000a; 2000b, who presented, in the 



same way, that disk-like systems are formed for a range 
of A = 0.10 — 0.04). However, the spin parameter of the 
virialized dark matter falls in a range 0.02-0.11 with a 
median value of 0.05, according to the results of N-body 
simulations (e.g., Barnes & Efstathiou 1987; Warren et al. 
1992) and analytical works (e.g., Steinmetz & Bartelmann 
1995; Eisenstein & Loeb 1995), in which A = 0.08 is near 
a high-end value. Kawata (1999) studied the evolution of 
a seed galaxy which has a slow rotation corresponding to 
A = 0.02 and showed that the end-product reproduces the 
observed properties of elliptical galaxies very well. Cos- 
mologically, a small spin-parameter is preferred in a galac- 
tic halo collapsing at a high redshift (Heavens & Peacock 
1988). Moreover, elliptical galaxies are considered to be a 
system whose stellar component is old and rotates slowly. 
Thus, it is a quite natural consideration that most elliptical 
galaxies were formed in a halo collapsing at high redshift 
and rotating slowly. 

Kawata's study was based on the hierarchical cluster- 
ing scenario, because in his simulation the galaxy forms 
through mergers of sub-clumps caused by initial small- 
scale density fluctuations (see Fig. 1 of Kawata 1999). 
The question we have to ask here is how the clustering 
of sub-clumps acts on the formation of elliptical galax- 
ies from a low-spin seed galaxy. In order to understand 
further the role of the clustering, we compare the evolu- 
tion of a low-spin seed galaxy including small-scale density 
perturbations with that of a uniform one, as the opposite 
case. Then, we focus on massive seed galaxies, and ex- 
amine whether both types of seed galaxy can reproduce 
the observed properties of BEGs quantitatively, because 
the properties of BEGs, as mentioned above, are expected 
to be able to provide constraints on their formation sce- 
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Fig. 1. — Time evolution of the system in Model 1. The upper, middle, and lower panels show the distributions of the dark matter, gas, 
and stars respectively. Each panel measures 150 kpc across and shows the x-z projection of the particles, where we set the z-axis to be the 
initial rotation axis. 



nario. If the uniform low-spin seed galaxy also succeeds 
in forming a system similar to BEGs, small-scale density 
fluctuations superposed on the initial slow-rotating sphere 
will not be indispensable to the model. The present study 
will thus clarify the role of the clustering of the sub-clumps 
induced by small-scale density fluctuations in galaxy for- 
mation from a massive low-spin seed galaxy. 

The plan for the remainder of the paper is as follows. 
In Section 2, the numerical method and the model are 
described. Section 3 presents results of numerical simula- 
tions. Conclusion is described in Section 4. 



2. METHODS 
2.1. The Code 

Details of our code are described in Kawata (1999). It 
is essentially based on the TrceSPH (Hernquist & Katz 
1989; Katz, Weinberg, & Hernquist 1996), which combines 
the tree algorithm (Barnes & Hut 1986) for the computa- 
tion of the gravitational forces with the smoothed par- 
ticle hydrodynamics (SPH; Lucy 1977; Gingold & Mon- 
aghan 1977) approach to numerical hydrodynamics. The 
dynamics of the dark matter and stars is calculated by 
the N-body scheme, and the gas component is modeled 
using the SPH. It is fully Lagrangian, three-dimensional, 
and highly adaptive in space and time owing to individual 
smoothing lengths and individual time steps. Moreover, it 
self-consistcntly includes almost all the important physi- 
cal processes in galaxy formation, such as self-gravity, hy- 
drodynamics, radiative cooling, star formation, supernova 
feedback, and metal enrichment. 

The radiative cooling which depends on the metallicity 
(Theis, Burkert, & Hensler 1992) is taken into account. 
The cooling rate of a gas with the solar metallicity is larger 
than that for a gas of the primordial composition by more 
than an order of magnitude. Thus, the cooling by metals 
should not be ignored in numerical simulations of galaxy 
formation (Kaellander & Hultman 1998; Kay et al. 2000). 

The star formation occurs where the gas density is 
greater than a critical density, p cr i t = 7 x 10 _26 g cm~ 3 , 



and the gas velocity field is convergent. The star forma- 
tion rate (SFR) is given by 



dp* 
dt 



dpg_ 
dt 



(i) 



where c» = 0.5 is a dimensionless SFR parameter and 
t g = y/3n/16Gp g is the dynamical time, which is longer 
than the cooling time in the region eligible to form stars. 
This formula corresponds to a Schmidt law in which SFR 
is proportional to p g - 5 . We apply a Salpeter IMF with 
the lower mass limit of 0.1M Q and the upper mass limit 
of 60M Q . After stars have formed, the massive stars 
(> 8Mq) die immediately, accompanied by supernova ex- 
plosions which inject thermal energy and metals to the sur- 
rounding gas. We adopt the same stellar nucleosynthesis 
yields as Kodama & Arimoto (1997). Type la supernovae 
are not taken into account in this paper. 

2.2. The Model 

We use the models which are almost the same as Kawata 
(1999). We consider an isolated sphere as a seed galaxy 
and compare the two models whose initial conditions are 
different. The only one difference between the two mod- 
els is whether small-scale density fluctuations following a 
CDM power spectrum are superimposed on the initial iso- 
lated sphere (hereafter Model 1) or not (hereafter Model 
2). In Model 1 the small-scale CDM density fluctua- 
tions are generated by Bertschinger's software COSMICS 
(Bcrtschinger 1995). To incorporate the effects of fluctu- 
ations with longer wavelengths, the density of the initial 
sphere has been enhanced and a rigid rotation correspond- 
ing to a spin parameter, A, has been added. The initial 
condition is determined by four parameters A, M to t, 08, in, 
and z c : the spin parameter is defined by 



A : 



J\E\ 



1/2 



GM, 



5/2 ' 



(2) 



where J is the total angular momentum of the system, 
E is the total energy, and M tot is the total mass of this 
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Fig. 2.— The same as Fig. 1 but for Model 2. 



sphere, which is initially composed of dark matter and 
gas; as, in is defined only in Model 1 and specifies an rms 
mass fluctuation in a sphere of radius 8 h~ x Mpc, which 
normalizes the amplitude of the CDM power spectrum; z c 
is the expected collapse redshift. If the top- hat density 
perturbation has an amplitude of Si at an initial redshift, 
Zi, we obtain z c = 0.36(5^(1 + z;) — 1 approximately (e.g., 
Padmanabhan 1993). Thus, when z c is given, Si at Zi is 
determined. Since we focus on the evolution of a massive 
low-spin seed galaxy, we set these parameters, A = 0.02, 
M to t = 4 x 10 12 Af Q , and z c — 2.07, which corresponds to 
an over-density of 2.5a for a CDM power spectrum with 
a s = 0.5 (White & Frenk 1991). In Model 1 we employ 
fs.in = 0.5. In Model 2 the initial sphere does not have 
any small-scale density fluctuations, i.e., the sphere has 
a completely uniform density. We assume a flat universe 
(fl = 1) with a baryon fraction of Clf, = 0.1 and a Hubble 
constant of Hq = 50 km s _1 Mpc -1 . 



angular momentum is not a specific one but the total value for each 
component. 



We carry out these two simulations using 9171 particles 
for gas and dark matter respectively. The masses of indi- 
vidual gas, stellar, and dark matter particles are 4.36 x 10 7 , 
4.36 x 10 7 , and 3.93 x 1O 8 M respectively, and softening 
lengths of gas, stellar, and dark matter particles are 2.06, 
2.06, and 4.28 kpc respectively. 
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Fig. 3. — Time variations of the star formation rate (upper panel), 
angular momentum (middle panel), and number of stellar clumps 
(bottom panel) in Model 1. In the middle panel, the thick solid line 
represents the evolution of the angular momentum of the baryon (gas 
and stars), and the solid and dotted lines show the time variation 
of the angular momentum of the gas and stars respectively. The 



3. RESULTS 

3.1. Time Evolution 

We simulated the evolution of the two models from 
Zi = 40 to z = 0. Figs. |l] and || show the evolution of 
the system in Model 1 and Model 2 respectively. These 
panels show the projection of the particles onto the x- 
z plane, where we take the z-axis to be the initial rota- 
tional axis. The evolution in Model 1 resembles Fig. 1 of 
Kawata (1999), although he simulated the evolution of a 
lower mass system (M tot = 8x 10 11 M ) than that of this 
paper. In Fig. 1 the system has already turned around at 
z — 3.21. We can observe the sub-clumps caused by initial 
small-scale density fluctuations. At z = 2.31, these clumps 
merge at the center of the system. The whole system has 
collapsed and settled into a spherical system at z = 1.89. 
After that, the system evolves little morphologically. The 
evolution of Model 2 is quite different (Fig. ||) . The gas is 
accreted to the center very smoothly, which induces star 
formation. Although the initial angular momentum is very 
small, a disk-like object is formed by contraction along the 
initial rotation axis. Model 2 also evolves little morpho- 
logically from z = 1.89 to z = 0. 
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Fig. 4. — Time variations of the star formation rate (upper panel) 
and angular momentum (lower panel) in Model 2. In the lower panel, 
the thick solid line represents the evolution of the angular momen- 
tum of the baryon (gas and stars), and the solid and dotted lines 
show the time variation of the angular momentum of the gas and 
stars respectively. The angular momentum is not a specific one but 
the total value for each component. 

Figs. H and [| show the time variation of the star forma- 
tion rate, the angular momentum of the baryon compo- 
nent, and the number of the stellar clumps in Model 1 and 
Model 2 respectively. Here, a clump means a bound ob- 
ject which is identified by applying to the stellar particles 
the friend-of-friend (FOF) algorithm with a linking length 
of 0.2 times the particle separation required to attain the 
mean density of the universe at each redshift and with a 
threshold particle number of 10. Since Model 2 initially 
has no internal small-scale density fluctuation causing for- 
mation of sub-clumps and its subsequent structure is also 
very smooth (Fig. g), we do not plot the history of the 
number of the clumps in Fig. ^. 

In Model 1 star formation begins in the sub-clumps be- 
fore the whole system collapses. The major star formation 
occurs due to the collapse of the whole system and contin- 
ues for about 1 Gyr. The star formation ceases owing to 
consumption of the gas. When the whole system collapses, 
the number of the clumps decreases suddenly (see also Fig. 
[j]) . This clearly shows that the galaxy is built through the 
mergers of sub-clumps, i.e., the hierarchical clustering. At 
the same time, the angular momentum of the gas com- 
ponent decreases well before the gas is transformed into 
stars. This provides an evidence that the mergers of the 
sub-clumps lead to the angular momentum transfer from 
the gas component to the dark matter. We have confirmed 
that the angular momentum of the dark matter increases, 
at the expense of the decrease in that of the baryon com- 
ponent, which consists of gas and stars. We can see that 
the formation of a spherical system supported by random 
motions in Model 1 is caused by the mergers of the sub- 
clumps which lead to the angular momentum loss of the 
baryon component and increase the velocity dispersion. 

On the other hand, in Model 2, the angular momen- 
tum of the baryon component is almost conserved during 
the collapse of the system (Fig. ||). The star formation 
begins just when the system has collapsed and hence the 
gas density has become high enough. For this reason, the 
star formation has a very high rate and a short period. 
Model 2 produces a disk-like object, because the baryon 
component conserves the angular momentum unlike Model 
I. Its collapse along the radial direction is halted at the 
point where the centrifugal force is balanced by the gravi- 



tational force, whereas the system is able to collapse along 
the initial rotation axis without limit. 

In the next section we present the properties of the end- 
products of the above evolution. In addition, we quanti- 
tatively examine whether each model can reproduce the 
observed properties of BEGs. 

3.2. Photometric and Kinematic Properties at z = 

Fig. H shows final distributions of stellar particles. These 
panels show a clear difference between the two configura- 
tions. The end-product of Model 1 has a nearly spherical 
shape, whereas that of Model 2 is flat. These stellar par- 
ticles contain the information about their age and metal- 
licity owing to the self-consistent calculation of the chemi- 
cal and dynamical evolution. By means of the population 
synthesis, we can derive the photometric properties of the 
stellar system from this information. Here, the spectral 
energy distribution (SED) of each stellar particle is as- 
sumed to be that of a single stellar population (SSP) that 
means a coeval and chemically homogeneous assembly of 
stars. The SED in a projected region, corresponding to 
one pixel of CCD, is obtained by the summation of the 
SEDs of the stellar particles whose projections fall within 
this region as follows: 



(3) 



where m„ U, and Zi are the mass, age, and metallicity of 
i-th stellar particle and f\(t,Z) is the SED of SSP of I 
M Q with the age, t, and the metallicity, Z. Here we used 
the data of SSPs of Kodama & Arimoto (1997). Since 
the observational data with which our results should be 
compared provide the luminosity distribution projected to 
a plane, we first derive the projected distribution of SED 
from the three dimensional distribution of stellar parti- 
cles. Fig. U shows the B band images corresponding to 
the projected distribution of stellar particles. Here a 501 
x 501 pixel mesh is chosen to span the squared region 
with 50 kpc on a side, and the flux of each stellar particle 
is smoothed using a gaussian filter with the filter scale, rps, 
of 1/4 times the softening length of the stellar particle, i.e., 
»*fs = 0.25e s . We discuss the dependence of the results on 
rps later. These images provide quite similar information 
to the imaging data obtained in actual observations. Thus 
we can obtain various photometric properties from these 
images in the same way as in the analysis of observational 
imaging data. 

Fig. ^ shows the azimuthally averaged surface- 
brightness profiles obtained from the x-y projected B band 
image shown in Fig. |^. In getting these profiles, we set the 
center to the position of a pixel which has the maximum 
V band luminosity. We set the width of each annulus (bin 
size) which defines the sampling range to half the soft- 
ening length. The error bars correspond to the standard 
deviation within each annulus. Since the error is calcu- 
lated in the flux, like <Tj b = — (fs) 2 , the error in the 
magnitude is defined, using the Taylor expansion, as 
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Fig. 5. — Final distribution of stellar particles (left) and B band image (right) for Model 1 and Model 2. The upper (lower) panels show the 
x-y (x-z) projections, where we take the z-axis to be the initial rotational axis. Each panel measures 50 kpc across. The contours correspond 
to the isophotes of 25, 23, 21, and 19 mag arcsec -2 respectively. 




Fig. 6. — The B band surface brightness profiles in Model 1 (left) and Model 2 (right). The solid lines denote the r 1 /™ law best-fit for the 
data plotted as the solid symbols. In the left (right) panel, the dotted line denotes the r 1 / 4 (exponential) law best-fit for the solid symbols. 
The fitting parameters are shown in the lower left corner of each panel. The error bars are shown only for the solid symbols and correspond 
to the standard deviation. 



Table 1 

Model final properties in the B band. 



Model 


log(Ie) 


n 




log(/ e ) or log(/ ) a 


r c or r D b 


M B 


V 




V/a 




(Lq,b PC 2 ) 




(kpc) 


(i QjB pc" 2 ) 


(kpc) 


(mag) 


(km/s) 


(km/s) 




1 


1.45 


4.66 


8.90 


1.51 


8.41 


-21.37 


20 


167 


0.12 


2 


2.64 


1.30 


3.08 


3.44 


1.89 


-21.38 


350 


217 


1.61 



a For Model 1, log(/ e ) of the r 1 / 4 law is listed. For Model 2, log(i"o) of the exponential law is listed. 
b For Model 1, r c of the r 1 / 4 law is listed. For Model 2, r^> of the exponential law is listed. 




_FlG. J. — Final B band images (left and middle) and the B band surface brightness profile (right) in Model 2. These are the same as Figs, 
pi and H, but the flux of each stellar particle is smoothed with a smaller rpg (for details, see text). 
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Fig. 8. — Final distribution of stellar particles (left), B band image (middle), and the B band surface brightness profile (right) in Model 2s. 
These are the same as Figs, a and pi but for Model 2s. 



These surface brightness profiles can well be fitted by the 
Sersic (r 1 /™) law, 

I(r) =I e \Q{-^[{r/r^'^l]}^ (5) 

where we adopt b n = 0.868n — 0.142, so that the effec- 
tive radius, r e , equals the half light radius in the range 
0.5 < n < 16.5 (Caon et al. 1993) and I e is the surface 
brightness at r c . This corresponds to the de Vaucouleurs 
(r 1 ' 4 ) law when n = 4, whereas the profile with n = 1 
corresponds to the exponential law, 

J(r) =I exp(-r/Rn), (6) 

where i?D is the scale length and Io is the surface bright- 
ness at the center. The solid lines in Fig. 2 show the 
best-fits in applying the r 1//n law. The dotted line in 
the left (right) panel of Fig. ^ shows the best-fit in ap- 
plying the de Vaucouleurs (exponential) law, i.e., holding 
n = 4 (n = 1) in equation (g). In fitting, we use the data 
with [is < 25.5 mag arcsec -2 . The best-fit parameters 
are shown in the lower left corner in each panel as well 
as m Table [TJ. The best-fit n in Model 1 is almost equal 
to four, and the de Vaucouleurs law also provides a good 
fit. Therefore, Model 1 succeeds in reproducing the light 
profile of BEGs, which is well fitted by the de Vaucouleurs 
law. On the other hand, the best-fit n in Model 2 is close 
to one, i.e., its profile is well described by the exponential 
law. This profile resembles the luminosity distribution of 
a disk galaxy rather than that of BEGs. However, the cen- 
tral luminosity is too large for disk galaxies (mo,B ~ 18.45 
mag arcsec -2 ). 

To see the robustness of these surface-brightness profiles, 
we examined the effects of rps and the number of particles 
on these profiles. Here we focus on Model 2, because Model 
2 provides a more compact system and is likely to be af- 
fected by the spatial resolution more severely than Model 
1 and the robustness of the profiles in Model 1 has already 
been examined in Kawata (1999). A large tfs is likely to 
lead to an artificially spread luminosity distribution and 
a large half-light radius. Fig. ^ shows the B band images 
obtained with rps = 0.125e s , i.e, half of rps in Figs. || and 
||. The radial surface-brightness profile obtained with a 
bin size of 0.25e s are indicated by points. Although the 
B band surface-brightness distribution is not smooth, the 
best-fitting parameters to the profile agree well with the 
ones in Fig. [| This demonstrates clearly that the profile 
in Fig. H is not artificially spread. In addition, we carried 



out the same simulation as Model 2 using a smaller num- 
ber of particles, i.e., 5575 particles for gas and dark matter 
respectively. We call this simulation "Model 2s" . The soft- 
ening lengths of gas, stellar, and dark matter particles are 
set to 1.82, 1.82, and 5.05 kpc respectively. Since Model 
2s employs a smaller number of particles and smaller soft- 
ening lengths for gas and stars, this model has a more 
favorable condition for the artificial two-body scattering 
than Model 2. Fig. ^ shows the luminosity distributions 
in Model 2s. Here we use tfs = 0.25e s . The artificial two- 
body scattering is expected to lead to a more extended 
surface-brightness profile than that of Model 2. However, 
there is little difference in the images and the profiles be- 
tween Model 2 and Model 2s. Therefore we conclude that 
the numerical resolution in Model 2 is sufficiently fine for 
the discussion of the luminosity distribution to be mean- 
ingful. 

It is known that a violent relaxation (Lyndcn-Bcll 1967) 
in sufficiently deep central potentials reproduces the de 
Vaucouleurs law (Hjorth & Madsen 1991 and references 
therein). A deep potential causes a strong scattering of 
tightly bound particles passing near the center, and leads 
to the outer envelope of the de Vaucouleurs profile. In 
addition, the dumpiness can strengthen the violent relax- 
ation so that particles are thrown to large radii, van Al- 
bada (1982) showed that a cold collapse with the formation 
and destruction of substructures is preferred in producing 
a de Vaucouleurs profile in dissipationless collapse. Carl- 
berg, Lake, & Norman (1986) found that dissipation, lead- 
ing to a deep central potential, could help get a good de 
Vaucouleurs profile. Katz & Gunn (1991) and Katz (1992) 
presented that dissipational hierarchical clustering is suit- 
able for bulge formation. Our results also indicate that a 
collapse accompanied with the dissipation and the mergers 
of the sub-clumps is quite appropriate for reproducing the 
de Vaucouleurs profile. 

Fig. shows the color and luminosity weighted metallic- 
ity gradients of the end-products in both the models. The 
color gradients are obtained by setting annuli of various 
radii in the B and R band images and subtracting the R 
band magnitude from the B band magnitude in each annu- 
lus. We made the projected images also for the metallicity, 
and obtained their radial profile. Points indicate the aver- 
age value in each annulus in the x-y projection. The width 
of each annulus is set to half the softening length. The 
error bars show the standard deviation in e ach annulus . 
The error in the color is written as cb-r = \ fL + cr 2 M ■ 



Fig. 9. — Color (upper panel) and metallicity (lower panel) gradients in Model 1 (left) and Model 2 (right). The solid lines show the best-fit 
linear relations for the data plotted as the solid symbols. The gradients in the best-fit lines are shown in the lower left corner of each panel. 
The error bars are shown only for the solid symbols and correspond to the standard deviation. 



Table 2 

Color and Metallicity Gradients. 



Model A(B — R)/A log(r) Alog(Z/Z©)/Alog(r) 



1 -0.127 -0.304 

2 -0.229 -0.776 



The profiles of color and metallicity [defined as log(Z/Z & )] 
are fitted by linear regression. In fitting, we excluded the 
data at radii less than the softening length and greater 
than the radius at which the B band surface brightness is 
\ib = 24.5 mag arcsec -2 (open symbols in Fig. ||), because 
the inner region is affected by the smoothing of gravita- 
tional forces whereas the number of particles within an 
annulus is too small in the outer region. The best-fit gra- 
dients are shown in the lower left corner of each panel in 
Fig. ^ and Table |2|. The data obtained from the simulation 
results have large errors leading to large errors in the cal- 
culated gradients. However negative gradients are clearly 
seen in colors and metallicity, i.e., the color (metallicity) at 
the center is redder (higher) than that in the outer region. 
From the observational studies, it is known that typical 
BEGs have A(B - i?)/Alog(r) = -0.09 ± 0.02 (Peletier 
et al. 1990) and Alog(Z/Z Q )/Alog(r) = -0.30 ± 0.12 
(Kobayashi & Arimoto 2000 and references therein). In 
Model 1 the gradients of both the color and the metallicity 
are in fairly good agreement with the observed gradients 
of BEGs. On the contrary, these gradients are too steep in 
Model 2. Model 1 is thus favored over Model 2 in repro- 
ducing the color and metallicity gradients as well as the 
surface-brightness profiles. 
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Fig. 10. — The half-light radius (upper panel) and V/cr (lower 
panel) against the absolute magnitude of the galaxy. The results of 
Model 1 and Model 2 are plotted as a circle and a square respec- 
tively. Crosses indicate the data for the ellipticals in the Virgo and 
Coma clusters (Faber et al. 1989; Davies et al. 1983). 

Carlbcrg (1984) showed that the dissipational collapse 
of a proto-galaxy with a uniform density leads to a steep 
metallicity gradient of Alog(Z/Z Q )/Alog(r) ~ —0.5, us- 
ing the numerical simulation. It has already been a se- 
rious problem that a metallicity gradient expected in the 
uniform collapse becomes steeper than the observed one 
(e.g., Davies, Sadler, & Peletier 1993; Kobayashi & Ari- 
moto 2000). The negative gradients in both the color and 
the metallicity mean that the star formation persists over 
longer time in the inner region than in the outer region, 
due to dissipative infall of gas. Since the residual gas is 
polluted by the past star formation, metal-rich stars are 
formed in the central region (Kawata 1999). In Model 1 
the star formation which occurs in the sub-clumps before 
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the whole system collapses decreases the amounts of the 
residual gas infalling dissipatively, leading to an appropri- 
ate shallower gradients. In addtion, the stellar mergers 
are expected to weaken an existing gradient (White 1980). 
Davies et al. (1993) supposed that BEGs have to experi- 
ence some stellar merger events. The scatteringinduced by 
the clustering of the stellar sub-clumps (Fig. |l|) in Model 

1 is also likely to play a quite important role of diluting 
the metallicity gradient which leads to the color gradient. 
The shallower gradients of Model 1 seem to be caused by 
the combination of these two effects. The results of our 
realistic numerical simulations suggest that the formation 
and clustering of sub-clumps are preferred in producing 
the observed color and metallicity gradients. 

In Fig. [h] we compare the sizes and the kinematical 
properties of the final stellar systems with those of the 
observed BEGs. The upper and lower panels respectively 
plot the half-light radii, r 0j s, and the ratios of the rotation 
velocity to the average velocity dispersion within 0.5r Ct s, 
V/a, against the B band absolute magnitude. The data of 
elliptical galaxies in the Virgo and Coma clusters are plot- 
ted as crosses, where we used the total magnitudes and the 
half-light radii of Faber et al. (1989), and V/a of Davies et 
al. (1983), assuming distance moduli of 31.0 for the Virgo 
(Graham et al. 1999) and 34.7 for the Coma (Bower et al. 
1992). There is no difference between Model 1 and Model 

2 in the total magnitudes of the final objects. However, 
the half-light radius in Model 1 is much larger than that 
in Model 2 and is in good agreement with that of the ob- 
served elliptical galaxies with the same luminosity. As to 
the kinematic properties, Model 1 succeeds in reproducing 
a low V/ a of BEGs very well. Owing to the loss of the an- 
gular momentum (Fig. ||), the end-product in Model 1 is 
a system supported by the random motions and hence ro- 
tating slowly. Furthermore, the small-scale initial density 
fluctuations allow star formation in the sub-clumps before 
the whole system begins to collapse and the collisionless 
property of the stellar component prevents the dissipative 
contraction of the system, leading to an appropriate half- 
light radius in Model 1. On the other hand, in Model 
2 the end-product is a nearly rotation-supported system, 
owing to the angular momentum conservation (Fig. 0). 



In addition, the initial-small angular momentum requires 
a large contraction along the radial direction for produc- 
ing enough centrifugal force, leading to a small half-light 
radius. These results demonstrate that, while Model 2 
is clearly unacceptable, Model 1 is pretty appropriate for 
producing BEGs. 

4. CONCLUSION 

We showed that the luminosity profile, color and metal- 
licity gradients, half-light radius, and V/a of the end- 
product formed by clustering of sub-clumps are in very 
good agreement with those of the observed BEGs quan- 
titatively. Since the CDM cosmology naturally predicts 
the existence of massive low-spin seed galaxies including 
small-scale density fluctuations, the CDM scenario is quite 
suitable for the formation of BEGs from a massive low-spin 
seed galaxy, in which the clustering of sub-clumps plays an 
essential role. 

On the other hand, in a massive low-spin seed galaxy in- 
cluding no initial small-scale density fluctuation, the small 
initial spin leads to a compact disk-like object dominated 
by the rotation owing to the angular momentum conser- 
vation. The disk has an exponential light profile, which 
is obviously inconsistent with the light profiles of the ob- 
served BEGs. Furthermore, the central surface brightness 
is too high to be in agreement with that of disk galaxies. 
No object like this has been observed in our universe de- 
spite its high surface brightness. This result suggest that 
a massive low-spin seed galaxy is not allowed to have a 
uniform density. 
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